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ABSTRACT  This  work  improves  Denton's  finite  volume 
time  marching  methods.  It  takes  this  type  of  method,  in 
which  there  is  1st  order  accuracy  in  time  and  2d  order 
accuracy  in  space,  and  improves  it  to  be  a  calculation 
method  in  which  there  is  2d  degree  accuracy  in  time  and  2d 
degree  accuracy  in  space.  In  conjunction  with  this, 
application  is  made  to  numerical  simulations  of  inlet  flow 
fields  associated  with  the  interior  flows  of  S  type  inlets 
as  well  as  exterior  flows  which  inlets  possess. 
Calculations  are  carried  out  for  actual  cases,  and  we 
obtained  reasonable  calculation  results. 

KEY  WORDS  Inlet,  Numerical  Calculation,  Numerical 
Simulation 


I .  FORWARD 

Using  time  marching  methods  to  solve  original  parameter 
equations  is  an  effective  method  to  make  stable  solutions 
for  non-viscous,  rotational  flow  fields.  It  is  possible  to 
very  easily  solve  subsonic,  transsonic,  and  supersonic  flow 
fields.  As  far  as  the  use  of  Euler  equations  to  solve  flow 
fields  is  concerned,  it  is  only  possible  to  carry  out 
nonviscous  numerical  simulations.  However,  it  is  possible 
to  make  preliminary  analyses  of  flow  fields  to  a  reasonable 
level.  At  the  same  time,  Euler  equation  solutions  are  also 
the  foundation  for  numerical  value  calculations  involved 
with  viscosity.  Denton's  methods ' are  a  type  of  time 
marching  method  associated  with  finite  volumes  and  implicit 
equations  of  2d  order  accuracy.  They  are  primarily  applied 
to  flow  field  calculations  for  blade  devices,  and  are  a  type 
of  effective  and  reliable  calculation  form.  The  work  in 
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question  takes  this  type  of  method  and,  from  1st  order 
accuracy  in  time  and  2d  order  accuracy  in  space,  improves  it 
to  a  calculation  method  with  2d  degree  accuracy  in  time  and 
2d  degree  accuracy  in  space.  In  conjunction  with  this,  it 
is  applied  to  numerical  simulations  associated  with  flow 
fields  of  S  inlet  internal  flows  as  well  as  incoming 
internal  and  external  flows.  It  achieves  reasonable 
results . 


II.  EULER  EQUATIONS  WHICH  CONSERVE  FORMS 


The  Euler  equation  set  which  conserves  forms  is 


dU 

dt 


+ 


dF  dG 
dx  dy 


(1) 


In  this  A170 

U=^{p,  PF„  pFy,  P-J)' 

F=(pv..  pi'iFp,  py.y„  py.h')' 

G^(pyy,  pyj'y,  pyi+p,  py^n' 


Corresponding  status  equations  are. 

(2) 

h*==c,r+-j-(yji+y}) 

e  =h*-Plp  (3) 

P  =PRT 


(4) 
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Talcing  Euler  equation  (1)  and  applying  it  to 
infinitesimal  control  volume  V,  one  obtains  the  change 
^  TJ  through  time  increments  ^t 


(5) 


^  n 

Here  At/=  CAP,  A  (PFJ,  A(Pf^y),  A  (Pe)y 


From  Euler  equation  (1),  one  solves  steady  state  flow 
fields.  With  regard  to  uniform  incoming  flows,  it  is 
possible  to  assume  that  enthalpy  h*  is  a  constant, 
reducing  by  one  the  required  iterative  equations  (energy 
equations ) . 


III.  DENTON  CALCULATION  METHODS ^ ^ ^ ^ 


As  far  as  the  use  of  central  difference  forms  to  solve 
non-stable  state  Euler  equations  is  concerned,  it  is 
necessary  to  opt  for  the  use  of  smoothing  measures  in  order 
to  guarantee  the  stability  of  solutions.  The  instability 
characteristics  of  center  difference  forms  are  clearly  shown 
in  wave  motion  changes  associated  with  flow  field  parameter 
highs  and  lows  which  draw  the  nodal  points  on  flow  lines, 
that  is,  high  speed,  low  pressure  nodal  points  are 
consecutive.  In  order  to  resolve  these  problems,  it  is 
possible  to  opt  for  the  use  of  types  of  head  wind  forms.  When 
handling  momentum  equations,  infinitesimal  units  of  change 
AU  (Fig.l)  are  directly  added  to  nodal  points  down  the 
flow.  Moreover,  one  makes  the  pressures  associated  with 
nodal  points  down  the  flow  be  directly  applied  to  nodal 
points  up  the  flow,  that  is 
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Fig.l  Infinitesimal  Body 

+  -2— C  A  (Pf^r);.;+  A 

This  form  strengthens  the  relationships  between  flow 
field  parameters  for  nodal  points  down  the  flow  and  nodal 
points  up  the  flow.  By  pressure  restraints  associated  with 
nodal  points  down  the  flow,  the  speeds  of  nodal  points  up 
the  flow  are  moved  up  or  down.  The  accuracy  is  only  1st 
order.  It  is  possible,  with  regard  to  the  addition  of 
corrective  factor  CFP  values  to  pressure  guantities,  to 
raise  the  accuracy  of  forms. 
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(9) 


pr.)-pi.!+cFPrj 

CFPr.)^(l  -RF)CFP';,,  +  RF(pl,,r 


-Pl.i) 


(10) 


RF  is  a  relaxation  factor  (normally,  approximately  equal  to 
0.05).  At  times  when  one  reaches  stable  reductions  or  convergence 
restraint,  the  equations  above  become 


As  far  as  the  forms  above  going  through  corrections  in 
pressure  quantities  is  concerned,  this  makes  form  accuracy 
become  2d  order.  Moreover,  stability  characteristics  do  not 
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change.  This  form  is  called  Form  A. 

Due  to  the  nature  of  flow  fields,  at  times  with  low 
Mach  numbers,  changes  in  pressures  and  densities  are 
basically  consistent.  Because  of  this,  the  pressure 
corrections  in  Form  A  are  capable  of  changing  into  ones 
which  are  realized  from  mass  equation  forms,  that  is  , 
density  changes  in  infinitesimal  bodies  A  are  applied  to 
upper  flow  surfaces.  One  gets  Form  B. 


Pr =  PI. 1-1  +  +  AP;.,./) 


(12) 


Others  are  the  same  as  Form  A. 

Form  B  is  capable  of  being  used  in  low  Mach  number 
calculations  existing  in  back  flows  or  stagnant  points. 
However,  as  far  as  those  existing  in  high  Mach  number  areas 
of  flow  fields  are  concerned.  Form  B  will  be  unstable. 

With  regard  to  the  combining  of  Form  A  and  Form  B, 
it  is  possible  to  present  a  type  of  form  (Form  C) 
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(13) 


p7,) = ink .  ^cFRor.))  ■  RTr.) 

CFRO]:]=^{  I -RF)-  CFROli  +  RFiPli-pi,;,,)  ^  ^ 

;  =  (p^^.) -.,  + -|-C  A  (PJ^,)1.;+ A  (pf^,) 

(pf'",)  ■:/ = (py>)i,+-~CA(py,)'}.r+  A  (pf^y)i,.,:i 

(16) 


In  this,  CFRO  is  a  density  correction  factor. 

When  one  reaches  a  stable  solution,  equation  (13) 
becomes 

P‘.>^Pi.,RT,,, 


This  form  captures  shock  wave  energy  strengths.  In 
conjunction  with  that,  it  is  stable  with  regard  to  all  Mach 
numbers.  However,  it  is  not  capable  of  being  applied  to 
calculations  for  flow  fields  with  reflux  or  back  flows. 

Finally,  it  is  possible  to  take  Form  B  and  Form  C  and 
combine  them  into  a  unified  calculation  form,  that  is,  take 
equation  (12)  and  equation  (13)  and  change  them  into  the 
three  equations  below 


Pr)-.  =  P;,/.+ 


a  (AP1.,  +  AP?-,.;) 


Pr'  =  P?,/+-4“(  1  ~  a  )(AP7.;  +  AP^-..,) 


(17) 

(18) 


(19) 


In  the  equations,  a  is  a  distribution  function. 
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When  Ma  <  1 ,  a  =  1 
Ma  ^  1 f  a  =  0 

However,  giving  a  in  this  way,  when  solving  two  and 
three  dimension  flow  fields,  one  will  see  the  appearance  of 
continuity  problems.  a  is  capable  of  being  defined  as  the 
temperature  function 


a  >  0 


(20) 


*  t 

In  the  eguation,  is  the  overall  temperature.  T  is 
critical  temperature. 

Energy  eguations  are  capable  of  using  center  difference 
forms  or  processing  eguations  of  the  same  types  as  eguation 
(17)  and  equation  (18)  with  regard  to  continuity  equations. 

This  form  is  capable  of  calculations  in  various  types 
of  flow  fields  with  high  and  low  Mach  numbers  as  well  as 
those  with  reflux  or  back  flows. 

Distribution  function  a  does  not  influence  stable 
solutions  after  convergence.  It  only  guarantees  form 
stability  with  regard  to  flow  field  calculations.  When 
calculating  shock  waves,  one  must  still  introduce  2d  order 
artificial  viscosity  terms  in  order  to  eliminate  shock  wave 
oscillations.  Then,  equation  (19)  becomes 


+  (p;:/w  -  p;:')  (pvt).,  -  p-r/.j/P-D 


(21) 


Denton  forms,  in  reality,  are  a  type  of  mixed 
difference  form.  In  drawing  up  flow  line  directions,  with 
regard  to  Euler  equations,  one  should  use  reflux  difference 
forms.  In  the  case  of  all  vertical  directions,  one  should 
use  center  difference  forms.  The  results  are  1st  order 
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accuracy  in  time  and  2d  order  accuracy  in  space.  Below, 
corrections  are  carried  out  on  solutions  associated  with 
U.  It  is  possible  to  make  Denton  calculation  forms  rise 
from  1st  order  accuracy  in  time  to  2d  order  accuracy  in 
time . 

From  equation  (1),  one  obtains 

dU  _  /  dF  dG  \ 

at  ~  \  ax  ay  ) 


Now  make 


Because 


hU 


fl+1  _ 


aU  Y 
at  J 


At  + 


aF  ,  aG 
ay 


ax 


+  - 


AH- 


1  faF  ^  aGY 

2  \  dy  J 


At+  0  (A/’) 


equations  above  then  take  their  integrals  versus 
infinitesimal  bodies  A  V,  and  one  obtains 


=  - \-AU+  0  (A/’) 


With  equation  (22),  one  corrects  ,  A(  P  ,  A(p 
Vy),  and  £s{pe).  One  is  then  able,  presupposing  virtually 
no  increase  in  the  amount  of  calculations,  to  raise  the 
accuracy  of  the  forms. 

When  applying  the  equations  in  question  to  carry  out 
numerical  simulations,  iterative  processes  associated  with 
equations  should  be 
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(1)  In  solving  continuity  equations,  to  solve  for  new 
densities 

(2)  In  solving  energy  equations,  to  obtain  new 
internal  energies 

(3)  Using  new  densities,  new  internal  energies  as  well 
as  old  velocities,  to  calculate  new  pressures 

(4)  Using  new  pressures,  solving  momentum  equations  to 
get  new  speeds. 

Using  this  sequence  to  treat  Euler  equations,  one  is 
able  to  guarantee  optimum  stability  in  association  with  the 
forms . 


IV.  STABILITY 


Denton  methods  are  time  marching  methods  associated 
with  implicit  equations.  Their  time  increment  t  is 
limited  by  stability  conditions: 

+  a 


1  is  the  grid 
calculations,  it 
simplify  it  to  be 


infinitesimal  flow 
is  possible  to  take 
Ax 


Ai<FT 


^/KRT, 


direction  width.  In  actual 
the  above  equation  and 


In  this,  FT  is  the  time  increment.  Generally,  FT  =  0.1-'^’ 

0.5. 

In  actual  flow  field  calculations,  one  will  also 
normally  see  the  appearance,  in  iterative  processes,  of 
excessively  large  flow  field  parameter  changes  associated 
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with  certian  special  points,  producing  non-stable  state, 
unstable  phemomena .  As  far  as  the  appearance  of  this  type 
of  situation  is  concerned,  if  control  is  not  added,  it  is 
possible  that  it  will  lead  to  iterative  divergence.  One 
effective  method  for  resolving  this  problem  is  to  introduce 
into  the  procedure  feedback  factors.  They  are  only  applied 
when  there  is  a  trend  toward  the  appearance  of  instability. 
Reducing  the  corresponding  grid  point  parameter  changes 
restrains  the  development  of  unstable  influences, 
guaranteeing  iterative  convergence 

V.  BOUNDARY  CONDITIONS 

As  is  shown  in  Fig. 2,  within  S  inlets,  there  are  three 
types  of  boundaries  all  together:  intake,  exhaust,  and  wall 
surface  or  skin  boundaries.  In  Fig. 3,  what  is  shown  is  S 
inlets  possessing  incoming  flow  exterior  fields.  All 
together,  there  are  five  types  of  boundaries:  incoming 
flow,  distant  flow,  exterior  flow  exhaust  or  exit,  interior 
flow  exhaust  or  exit,  as  well  as  wall  surface  or  skin 
boundaries . 

It  is  possible  to  know,  from  deductions  of 
characteristics  theory,  that,  when  incoming  flows  are 
subsonic,  incoming  flow  boundaries  or  interior  flow  intake 
boundaries  should  fix  three  flow  movement  parameters. 
Generally,  they  fix  overall  temperature  T^,  overall 
pressures  p^ ,  as  well  as  flow  movement  angle.  When 
interior  flow  exhausts  or  exits  or  exterior  flow  exhausts  or 
exits  have  subsonic  flows,  their  corresponding  boundaries 
should  fix  one  flow  movement  parameter.  Generally,  they  fix 
static  pressure  distribution.  With  regard  to  wall  surface 
or  skin  boundaries,  it  is  possible  to  use  mutually  tangent 
flow  movement  conditons.  As  far  as  the  exterior  flow 
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distant  field  boundaries  of  Fig. 3  are  concerned,  it  is  only 
necessary  for  boundaries  to  be  sufficiently  distant  from 
inlets,  and  it  is  possible  to  fix  or  designate  free  flow 
autonomous  conditions. 


Fig. 2  S  Inlet  Interior  Flow  Grid  (1)  Intake  (2)  Upper 
Wall  (3)  Lower  Wall  (4)  Center  Line  (5)  Exhaust  or 
Exit 
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Fig. 2  S  Inlet  Grid  Possessing  Exterior  Flows  (1)  Incoming 
Flow  (2)  Distant  Field  (3)  Exterior  Flow  Exhaust  or  Exit 
(4)  Upper  Wall  Surface  or  Skin  (5)  Lower  Wall  Surface  or 
Skin  (6)  Center  Line  (7)  Interior  Flow  Exhaust  or  Exit 
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VI.  INITIAL  FLOW  FIELDS  AND  GRIDS 


What  is  capable  of  defining  initial  flow  fields  is: 

(1)  overall  temperature  and  overall  pressure  along  the  drawn 
up  flow  lines  being  constant  (2)  static  pressure  along 
drawn  up  flow  lines  changing  in  a  linear  way,  and  (3)  the 
flow  direction  of  the  entire  flow  field  conforming  to  the 
direction  of  a  drawn  up  central  flow  line.  Numerical  tests 
demonstrate  that  defining  initial  flow  fields  in  this  way  is 
feasible . 

As  shown  in  Fig. 2  and  Fig. 3,  grids  are  composed  of 
drawn  up  flow  lines  and  vertical  transversals. 

VII.  CALCULATED  INSTANCES  AND  CONCLUSIONS 

Fig. 4(a)  is  a  calculation  grid  for  an  interior  convex 

channel.  It  is  composed  from  31x13  individual  nodal  points. 

Fig. 4(b)  is  a  subsonic  calculation  result  for  a  channel  of 

this  geometry.  Flow  field  intake  Mach  numbers  are  all  0.5. 

Symmetry  characteristics  for  these  results  are  good.  They 

are  basically  the  same  as  finite  volume  time  marching  method 

( 3 ) 

calculation  results  associated  with  ^  ^^^^^^ear  )  ’ 

Fig. 4  (c)  is  the  convergence  history  for  this  calculated 

instance.  There  were  all  together  1000  iterative 

increments.  The  maximum  residual  difference  associated  with 

—  5 

X  direction  velocities  was  smaller  than  1.1x10 
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Fig. 4(a)  Convex  Channel  Calculation  Grid  (1)  Upper  Wall 
(2)  Lower  Wall 
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Fig. 4(b)  Mach  Number  Distributions  Associated  with  Channel 
Upper  and  Lower  Walls  (1)  Upper  Wall  (2)  Lower  Wall 


Fig. 4(c)  Convergence  History  (1)  Time  Increment  Number 


Fig. 5  Mach  Number  Distributions  Associated  with  Two 
Dimensional  S  Type  Interior  Flow  Channels  (1)  Center  Line 

(2)  Upper  Wall  (3)  Lower  Wall 
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Fig. 5  is  calculation  results  for  a  two  dimensional  S 
type  channel  as  shown  in  Fig. 2.  The  grid  points  are  23x7. 
Boundary  conditions  for  incoming  flows  are  given  as  uniform 
overall  pressure  p^  equal  to  101330.0  Pa,  uniform  overall 
temperature  of  288  K,  as  well  as  y  direction  velocity 
component  V  equal  to  zero. 

Fig. 6  is  the  results  gotten  from  carrying  out 
calculations  after  the  equivalent  section  associated 
with  a  10%  addition  in  length  to  the  Fig. 2  inlet  position. 
The  boundary  conditions  given  are  the  same  as  those  for 
Fig. 5  calculated  instances.  This  result  eliminates  the 
pulsations  in  flow  field  distributions  at  intake  positions 
in  Fig. 5  calculated  instances.  The  reason  is  that  one  is 
given  an  incoming  flow  boundary  condition  of  Vy.  =  0. 
Addition  of  equivalent  sections  makes  it  possible  to 
guarantee  that  the  realization  of  these  conditions  will  be 
even  more  rational. 


Fig. 6  S  Inlet  Mach  Number  Distribution  Associated  with  the 
Addition  of  Equivalent  Sections  (Not  Including  Exterior 
Flows)  (1)  Upper  Wall  (2)  Lower  Wall  (3)  Center  Line 
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Fig. 7  Mach  Number  Distributions  Associated  with  S  Type 
Inlets  (1)  Center  Line  (2)  Upper  Wall  (3)  Lower  Wall 

Fig. 7  is  calculation  results  for  S  type  inlets 
possessing  incoming  flow  exterior  fields  as  shown  in  Fig. 3. 
Incoming  flow  and  interior  flow  exit  boundary  conditions  are 
the  same  as  those  for  Fig. 5  calculated  instances.  In 
conjunction  with  thiS/  the  given  distant  field  boundary 
conditions  and  incoming  flow  boundary  conitions  are  the 
same.  Exterior  flow  exhaust  or  exit  boundary  conditions  are 
given  as  uniform  static  pressure  p^  egual  to  83688.4  Pa. 
Results  clearly  show  that  Fig.7's  Mach  number  distribution 
and  FigiG's  Mach  number  distribution  are  basically 
consistent.  This  type  of  result  is  rational. 

Going  through  several  calculated  instances,  one  comes  to 
a  number  of  conclusions. 

(1)  Denton  method  calculation  efficiencies  on  IBM  4341 
computers  are  each  grid  point  each  time  increment 
approximately  1.0xl0“^s.  One  normally  reguires  500 
iterative  increments.  Maximum  residual  errors  associated 
with  X  direction  velocity  components  in  flow  fields  are  all 

_3 

capable  of  being  less  than  0.5x10 
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(2)  Convergence  standards  can  be  fixed  to  be  maximum 

—  ? 

residual  errors  smaller  than  0.5x10 

(3)  Euler  equation  processing  procedures  influence 
solution  stability  characteristics. 

(4)  Denton  method  numerical  simulations  with  regard  to 
S  inlet  flow  fields  are  feasible.  Using  Euler  equations  to 
simulate  non-viscous  flow  fields,  it  is  possible  to  provide 
flow  field  distribution  reference  values. 
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C535  AVRAD0C9V1SARCCM 

1 

C539  TRASANA 

1 

Q592  ESTC 

4 

Q619  MSIC  REDSTONE 

1 

Q008  Nnc 

1 

Q043  AEyEC-IS 

1 

E051  HQ  USAF/INET 

1 

E404  AEDC/DOF 

1 

E408  AIWL 

1 

E410  AFDTC/IN 

1 

E429  SD/IND 

1 

POOS  DOE/ISA/DDI 

1 

P050  CIVOCEVADD/SD 

2 

1051  AFIT/IDE 

1 

P090  NSA/CDB 

1 

2206  ESL 

1 

Microfiche  Nbr:  FTD94C000492 
NAIC-ID (RS) T-0922-92 


